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In many cases of interest, the perturbative series based on conventional Feynman diagrams have 
a zero radius of convergence. Series with a finite radius of convergence can be obtained by either 
introducing a large field cutoff or by replacing the exponential of the perturbation by a sequence 
of approximants as recently proposed by Pollet, Prokof'ev, and Svistunov. We compare these 
two methods for integrals and quantum mechanical problems. The two methods perform well in 
complementary regime (strong coupling for the large field cutoff and intermediate coupling for 
the other method). We briefly discuss potential applications for lattice gauge theory with compact 
groups (which have a build-in large field cutoff). 
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1. Introduction 

More than 50 years ago, Dyson argued that the radius of convergence of the expansions of 
physical quantities for QED in powers of e 1 is zero because of the vacuum unstability if e 2 is 
changed to — e 1 [JTJ| _ Later, it was established that for many models, the perturbative expansions 
calculated from Feynman diagrams are divergent For instance, for several models with A0 4 
interactions, it is well-known that the perturbative coefficients increase at a factorial rate with the 
order. The failure of the convergence of the expansion can be understood easily from the simple 
integral where 

rv-^* 4 + £ t^. r d ^y^ . (i.D 

In this case, Dyson's instability is related to large field configurations [||, |J]. 

Several modifications of the large field configurations have been proposed. One possibility 
consists in truncating the large field contributions [gj: in the partition function of X<p 4 models, the 
integrals over (— oo,oo) are replaced by integrals over (—^ max ,^>max)- With this modification, Eq 
( |Ll| ) becomes an equality, but we are dealing with a different integral than the one we started with. 
On the other hand, for lattice gauge models with compact groups, there is a build-in large field 
cutoff and conventional series are obtained by "adding the tails of integration" 

Another modification consists in replacing the original partition function by a convergent se- 
quence of approximations where the exponential of the perturbation is replaced by a suitable sum 
of Bessel functions, as recently proposed by Pollet, Prokof 'ev, and Svistunov In the follow- 
ing, this method is abbreviated as the PPS method. As we will see, in this case, the contributions 
from the large fields have damped oscillation, and the series for the modified integrals become 
convergent. 

In this proceedings, we compare these two improved perturbative methods mostly for the 
simple integral A0 4 theory and mention recent progress for the anharmonic oscillator. 

2. Two Modified Perturbative Methods 

For models with A0 4 intereactions, and the perturbative expansion of the generating function 
of «-point functions can be written as a product of two exponents: 

Z[j] ~ e ~ kI * < '& 4 e^ jyG ™ A , (2.1) 

In the large field cutoff method, the second exponential is modified. In the PPS method, the first 
exponential is modified. 

2.1 The Large Field Cutoff Method 

When a large field cutoff is introduced, the standard shifts that allow to solve the free theory 
with sources need to be modified by writing the truncated integral as the integral over the whole 
real axis minus the integral over \<j>\ > ty max - We then obtain an expension of the form [§]: 

Z[j] ~ g - A ^37?) 4 e2^* G wA(l_c£ ( dO x e " ) - (; 1 ' + ...) . (2.2) 

y J ymax 
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If we want to write Feynman rules, the additional terms require a generalization of Wick theorem. 
Up to now, the modified coefficients have only been calculated numerically. 

2.2 The PPS Method 

In the PPS method, the vertex part is modified into a new sequence of approximants with a 
regulator N [0]. The partition function becomes 

Z[j]^Z N [j] = r d^ 2 [f{X^,N)] N 

J —DO 

* t X k C k , (2.3) 



k=0 

In which, /( X ? , N ) = J (a (N) X ? ) + 2J 2 {a (N) X ? ) + f J 4 (a (N) X ? ) , with a (iV) is chosen so 
that the series of f(X*(j>,N) has the form 1 - ^ + 0(^ft. 

By increasing the regulator iV, the modified partition function becomes more like the origi- 
nal one. Note that for the simple integral, \Ct\ is near-Gaussian and decays rapidly for large k. 
For a fixed X, the numerical value of the modified partition function is convergent as the order p 
increases. Feynman rules can be written easily: 

Z[j] - e -^^(^) 4+fc ( iV )(^^(35) 8 )+-] e |l y ^VG y ,J z . (2.4) 

The new vertices of order X 2 and higher can be seen as counterterms. For instance, the lowest order 
counterterm contributes to the vacuum energy at order X 2 through the diagram: 




3. Comparison of Modified Partition Functions in 0-D (simple integrals) 

In order to assess the quality of various perturbative expansions (or sequences of these), we 
ask the two following questions: 

• How well does the modified model approximate the original model? 

• How well does the perturbative series for the modified model at successive orders approxi- 
mate the modified model? 

These two questions can be addressed for a broad range of couplings and summarized in what we 
call "error graphs", namely the significant digits vs. logioX. In all the error graphs, the number of 
(correct) significant digits is calculated as 

modified value — true value 

-log io\ : . (3.1) 

true value 

In this section, we discuss these two questions for the simple integral (JTTTJ). We start with the first 
question. The error graph for the large field cutoff method with different cuts are shown in Fig.|I|. 
The error introduced by the field cutoff becomes smaller at strong coupling, because in this regime, 
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the large field contribution are strongly suppressed in the original integral. At small coupling, the 
error graph stabilizes at the error due to removing the tails in the Gaussian limit. In comparison, 
the error graphs for the PPS modified model with different regulator N are shown in Fig. |2| (Left). 
We calculated the asymptotic values for the PPS method at small and large A. At small A, the 
significant digits are — Zogio^xr- The asymptotic value for the PPS method at large A can be 
calculated numerically and are shown in Fig. || (right). The methods are complementary in the 
sense that they have chance to perform well in the opposite limits. 




Q I ' ' L_ I I I I 

-3 -2-10 1 



Log X 

Figure 1: Number of significant digits for the simple integral obtained from field cutoffs with (j) mclx = 2 
(blue), (j> max = 3 (red), mf „ = 4 (green). 

We can now address the second question regarding the approximation of the modified models 
by perturbative series. It seems clear that perturbative expansions of the modified models can at best 
do as well as the modified models (understood as solvable numerically with any desired accuracy). 
The results for the two modified methods are shown in Fig. ||. One sees that as the order increases, 
the error of the perturbative expansions gets closer to the error due to the modification. In other 
words, in both cases, the perturbative series converge toward the values taken for the modified 
models. Fig. || also includes regular perturbation theory and one can see that at sufficiently small 
coupling it provides the best answer. The field cutoff method works well at large coupling and the 
PPS method at intermediate coupling. 

4. Logarithms of the Modified Partition Functions in 0-D 

Unlike the perturbative series for the modified partition functions which are convergent for 
any complex value of A, the expansion of their logarithm have a finite radius of convergence which 
is the absolute value of the zeros of the modified partitions function closest to the origin in the 
complex A plane. However, accurate values can be obtained by using Pade approximants for both 
methods. As the order of the Pa.de approximant increases, the significant digits of the expansion 
can be improved and we reach the numerical limit corresponding to the modified models. 
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Figure 2: Left: Number of significant digits obtained from the PPS modification for the simple integral with 
regulator jV = 2 (blue), N = 3 (red), N = 4 (green). Right: Asymptotic value (black) for the PPS method at 
large X: N = 2 (blue), N = 3 (red), N = 4 (green). 




Figure 3: Number of significant digits for the simple integral obtained with a field cutoff with <j) max = 2 
at order 7, 11, 15, as a function of X (red lines with the blue asymptotic line increasing as X increases), 
from the PPS method with = 2 at the same orders (red lines with the blue asymptotic line decreasing as X 
increases), and from regular perturbative method at the same orders (black). 

5. Preliminary Results for the Anharmonic Oscillator 

In the last two sections, we found that both the PPS method and the large field cutoff method 
work in some region of A . It is now time to ask if these methods can be used in more complicated 
models. In this section, both methods will be extended to the case of the anharmonic oscillator (ID 
X<p 4 ). The physical quantity we used to compare with is the ground state energy of the anharmonic 
oscillator. Following the idea from Ref. |7|], we can solve the numerical Schrodinger equation to 
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Figure 4: Error graphs for the anharmonic oscillator. Black: regular perturbation theory at order 15. Blue 
line (decreasing when X increases): the PPS method at regular N = 4. Green: the PPS method at order 
15. Blue line (increasing when A increases): cutoff method for a cut = 2. Red: cutoff method at order 
9,13,17,21,25. 



get the ground state energy for two modified methods and also their perturbative series. However, 
in the PPS case, the approximation of the potential by polynomials is more subtle because the 
new modified potential has logarithmic singularities. We are still working on this case, and the 
preliminary results are shown in Fig. || which shows patterns similar to Fig. ||. 



6. Preliminary Results on the Average Plaquette in U(l) Lattice Gauge Theory 

In order to extend the modified methods to U(l) lattice gauge theory, we need to check the 
easiest quantity, the average lxl plaquette. On the lattice, there are volume effects discussed in 
Ref. We reproduce these volume effects on our error graphs (Fig. ||), in which the numerical 
error is estimated through a weighted 20-seed data [gj. There is a gap between these two graphs due 
to the volume corrections. The sensitivity to the volume effects indicates that it might be possible to 



extract a "perturbative envelope" as in Ref. [ 10] and explain it in terms of large field contributions. 
This is still in progress. 



7. Conclusions 



By introducing a large field cutoff or replacing the exponential of the perturbation by a se- 
quence of approximants, we can get finite radius of convergence for the perturbative series of the 
logarithms of the simple integrals. The application of these two modified methods in more compli- 
cated models, like the anharmonic oscillator and U(l) lattice gauge theory is in progress. 



Acknowledgments 

We thank Alexei Bazavov and Daping Du for valuables discussions on U(l) lattice gauge 



6 



Comparison of Improved Perturb ative Methods 



Haiyuan Zou 



U(l) Plaquette V=oo 
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Figure 5: Comparison of the weak coupling expansion at different order with numerical errors without vol- 
ume correction (left) and with volume correction (right). The numerical error is estimated using a weighted 
20-seed data. 
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